KU-TP 056 



Beta Function and Asymptotic Safety in Three-dimensional 

Higher Derivative Gravity 



Nobuyoshi Ohta 1 

Department of Physics, Kinki University, Higashi-Osaka, Osaka 577-8502, Japan 

Abstract 

We study the quantum properties of the three-dimensional higher derivative gravity. In par- 
ticular, we calculate the running of the gravitational and cosmological constants. The flow of 
these couplings shows that there exist both Gaussian and nontrivial fixed points in the theory, 
thus confirming that the theory is asymptotically safe. It is shown that the new massive gravity 
or f(R) gravity in three dimensions do not correspond to the fixed point within the approxima- 
tion that the coefficients of the higher curvature terms are not subject to the flow. The fixed 
point value of the cosmological constant is found to be gauge-independent, positive and small. 
We also find that if we start with Einstein term with negative sign, the fixed point only exists 
when the coefficient of the Einstein term has positive sign. 
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1 Introduction 



One of the long standing important problems in theoretical physics is to understand quantum 
properties of gravity. The usual Einstein gravity is known to be non-renormalizable in four and 
higher dimensions. If we add Ricci and scalar curvature squared terms, the theory becomes 
renormalizable, but then the unitarity is lost [1]. So this higher derivative gravity does not 
appear to make sense as a physical theory. 

Recently a very interesting proposal has been made that the addition of such higher order 
terms to three-dimensional gravity makes the theory unitary if the coefficients are chosen ap- 
propriately [2] (with "wrong sign" Einstein term). The usual Einstein gravity does not have 
any propagating mode, but the addition of these terms introduces propagating massive gravi- 
ton around flat Minkowski and curved maximally symmetric spacetimes [anti-de Sitter (AdS) 
and de Sitter (dS) spacetimes]. Another theory of massive graviton with Lorentz-Chern-Simons 
(LCS) term has long been known as topologically massive theory [3], which breaks the parity. 
The above new theory is an ordinary parity-preserving theory and is called new massive gravity. 
This is very interesting in that we have really dynamical theory of gravity that is unitary even 
though higher derivative terms are included. Since then, various aspects of the theory have 
been investigated. Linearized excitations in the field equations were studied in [4]. Unitarity is 
proven for Minkowski spacetime in [5, 6, 7], whereas it is discussed in [8] for maximally symmet- 
ric spacetimes. A complete classification of the unitary theory for the most general action with 
arbitrary coefficients of all possible terms is given in [9]. The partial result of unitarity condi- 
tion on the flat Minkowski spacetime was known for the usual sign of the Einstein theory [10]. 
Unfortunately the new massive gravity turned out to be non-renormalizable though the general 
theory with arbitrary coefficients for the quadratic curvature terms is renormalizable [6, 11, 12]. 

However, even if the theory is not renormalizable, it is possible that a theory has a nice 
ultraviolet property if it has an ultraviolet fixed point. A non-renormalizable theory may be 
made effectively renormalizable by a rearrangement of the perturbation series or by addition of 
higher derivative terms. However terms of finite order in the perturbation series then contain 
what appear to be unphysical singularities. Such unphysical singularities may be almost certainly 
avoided if the couplings approach a fixed point in the ultraviolet energy. This property is known 
as asymptotic safety [13], and it seems that this is the only way to make sense of gravity theory 
known to date. In particular the asymptotic safety is a wider notion than the renormalizability. 
Any theory will always have a fixed point at the origin. If this is the only suitable fixed point 
with ultraviolet critical surface of nonzero dimensionality, then the asymptotic safety requires 
that the couplings lie on this surface. In order for the trajectory of the coupling constants to 
hit the origin in the high energy, all couplings with negative dimensionality (in powers of mass) 
should vanish. These are precisely the non-renormalizable interactions, so this sort of theory 
must be renormalizable in the usual sense. 

It should be very interesting to examine what quantum properties the above theory of higher 
derivative gravity has since it is dynamical theory of gravity and preserves parity as the ordinary 
gravity theories. The renormalization group (RG) properties of four-dimensional gravity are 
studied, for instance, in [14, 15, 16]. Those of the topological gravity in three dimensions have 
been studied in [17], but as far as we are aware, these have not been examined for general higher 
derivative theory in three dimensions related to the above new massive gravity. Discussions of 
four-dimensional higher curvature theories are given in Refs. [18, 19, 20, 21]. A simple RG flow 
argument is used in the case of Horava gravity in [22], and also 3D f(T) gravity is analyzed in 
[23]. In this paper we would like to examine if the three-dimensional theory is asymptotically 
safe, and the theory at the fixed point has any relation with the above new massive gravity. 
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In the next section, we briefly summarize the Wilsonian RG approach to gravitational theory. 
Here we have to introduce the cutoff to regularize the theory, and must be careful about what is 
meant by the cutoff since the definition of the cutoff in general uses a metric which is dynamical. 
We just follow the standard way to use background method. In sect. 3, we introduce the gauge 
fixing and the corresponding Faddeev-Popov (FP) terms, and derive the quadratic part of the 
action in the fluctuations. In sect. 4, we then evaluate the resulting functional traces and derive 
one-loop beta functions. In sect. 5, we discuss what flows and fixed points we obtain. In our 
analysis, we restrict ourselves to the flow in the gravitational and cosmological constants since it 
appears that the coefficients of the higher derivative terms can be set to constant, and then the 
analysis becomes simple. This analysis is done for the usual and opposite signs of the Einstein 
term. In both cases, we find that there are Gaussian fixed point at the origin and nontrivial 
fixed points as the ultraviolet attractive points. Even without the beta functions for the higher 
curvature terms, our analysis indicates that the new massive gravity or a special case of f(R) 
gravity in three dimensions do not correspond to the fixed point. We also find that if we start 
with Einstein term with negative sign, the fixed point only exists when the coefficient of the 
Einstein term has positive sign. The final section 6 is devoted to discussions and conclusion. 



2 Wilsonian method for renormalization group 

In the Wilsonian RG, we consider that the effective action describing physical phenomenon at a 
momentum scale k can be thought of as the result of integrating out all fluctuations of the fields 
with momenta larger than k. We can regard k as the lower limit of the functional integration 
and call it the infrared cutoff. The dependence of the effective action on k gives the Wilsonian 
RG flow. 

There are several ways of implementing this. We follow Ref. [17] and suppress the contribu- 
tion of the field modes with momenta lower than k by modifying the low momentum end of the 
propagator and leaving all the interactions unaffected. 

We start with a bare action S[(p] for some fields <j>, and add a suppression term ASk[<p] 
which is quadratic in the field. We choose some differential operator O whose eigenfunctions 
ip n , defined by Oip n = \ n (f n , are taken as a basis in the functional space we integrate over: 

<f>( x ) = ^2<t>n<Pn{x)- (2.1) 
n 

The additional term is written as 

AS fc [0 = \f dx4>{x)R k {0)4>{x) = ^4%R k (K)- (2.2) 

The kernel R k (z) is introduced as a cutoff, and is chosen to be a monotonically decreasing 
function in both z and k, namely Rk(z) — > for z 3> k and Rk{z) ^ for z <C k. Define a 
/c-dependent generating functional of the connected Green functions by 

e~ Wk ^ = J V(/)exp - S[<f>] - AS k [<f>] - J dxJcfy , (2.3) 

and a modified A;-dependent Legendre transform 

T k [4>] = W k [J] - I dxJ<\> - AS k [4>] • (2.4) 
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In the limit of k — > oo, this functional tends to the usual effective action r [</>]. 
We assume that admits a derivative expansion 

oo 

r*[^ffi] = EEW n) ^ ( 2 - 5 ) 

n=0 i 

(n) (n) 

where g\ are coupling constants and 0\ all possible operators constructed with the field <p 
and n derivatives compatible with the symmetries of the theory. The index i is used to label 
different operators with the same number of derivatives. At the one loop, is given by 

T^ = S + ^Tr\og(S^+R k ), (2.6) 

where denotes the second variation of the bare action. We then obtain 

The factor k^jjr goes to zero for z > k 2 . One can obtain the one-loop beta functions from this 
functional equation. 

We apply this method to higher derivative gravity in three dimensions. Here one has to be 
careful about what one means by cutoff because the definition of a cutoff generally makes use 
of a metric which in turn is to be treated as a dynamical field. We follow [14, 17] and use the 
background method, and effectively replace the dynamical metric g^ v by a spin-two field h^ u 
propagating in a fixed background g^ v . The background metric can be used to unambiguously 
distinguish what is meant by long and short distances. 



3 Quadratic expansion of the action 

We consider the action 



= Z J d 3 xj- 



-9 



aR - 2A + aR 2 + (3R 2 , (3.1) 



where Z = 1/(1QttG) is the three-dimensional gravitational constant, a,f3, [i and a(= ±1) are 
constants, and Ao is a cosmological constant. 

We consider the action up to second order around the background spacetime 

9^=9^ + h fM/ , (3.2) 

where we keep the formulae for arbitrary dimension D for the moment with future application 
to higher dimensions in mind. Now the background g„ v is chosen to be a maximally symmetric 
spacetime with the curvatures 

- R _ — R . _ _ _ _ . . . 

Riiv = jj9nv, Rnu P \ = _ j-j [9fi P 9uX ~ 9n\9v P ), (3.3) 



where R = ± P ^ 2 ^ with the + sign for de Sitter and — sign for anti-de Sitter spaces, and I is 

h = r v Kv (3.4) 



the radius. We define 
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Here and in what follows, bar indicates that the quantity stands for the background, the indices 
are raised and lowered by the background metric g, the covariant derivative V is constructed 
with the background metric, and the contraction is also understood by that. 
We parametrize the metric in general dimensions as 



1 



1 



with 



V A /i^ = 0, r v hl v = 0, V A 6 = 0. 
A straightforward calculation then yields the quadratic terms [9] 



(3.5) 



(3.6) 



D(D-i 



R 



P □- 



D(D - 1 



-R 



—R(Da + (3) + a 



(D-l)(D-2) 



4D 2 



fjAnfj - 2hAJ □ (□ + -j^—[R) v + hA(a + jj^-jR) h 



A 
+ 2L 



D - 1 



D-2 



(3.7) 



where we have defined 



4(D-l)a + Dp 2{D-4) 



A = 



D 



-□ - 



D-2 
A^A -— ftx 



D{D - 2) 



(Da + f3)R- a, 



2D 2 



-(Da + /3)R 2 . 



(3.8) 



These field redefinitions cancel the Jacobian introduced in the path integral by the parametriza- 
tion (3.5). 

Our next task is to introduce the gauge fixing and the corresponding FP ghost terms. The 
BRST transformation for the fields is found to be 

5b9 P v = -$X[g pu V p c p + g pp V u c p + V p g pu c p ], 

5 BC p = -6Xc p V p c p , 

$bc p = id\B p , 

S B B P = 0, (3.9) 

which is nilpotent. Here 5X is an anticommuting parameter. The gauge fixing term and the FP 
ghost terms are concisely written as 

Cgf+fp/Z = i6 B [c^n + c){ X >i ~ ^B p )]/5X 

= -B,(D + c) X » - «■„(□ + C ) + ^V„V + R^c" + ^B,{U + C )B», (3.10) 

where c is a constant, the indices are raised and lowered with the background metric, and 



Xu = V p h p u - ^-^V u h, 



(3.11) 
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with a and b being gauge parameters. The unusual factor (□ + c) is introduced in (3.10) in 
order to be able to diagonalize the kinetic terms of the gravitational fields [see Eq. (3.17) below]. 
Eliminating the auxiliary fields B^, we find the quadratic terms in the gauge fixing are 



(D - 1)(4 - D-bD) 
2D 2 



D - 1\2 



D 



(4 - D - bDf 



R* 



R 



D - 1 



R\ , 



hn\(n + c +—^h 



16D 2 

Similarly we find the corresponding FP terms are given by 

R\„ R\/3-b 



Cfp 



V"{n + c)(n + -)v lt -s(n + c+^ 



2 D + D R)S 



(3.12) 



(3.13) 



where a factor Z is dropped for simplicity (we can think that it is absorbed into the fields) and 
we have defined 



c M = ^ + v M s, s = V^s, v li v» = o, 
c^ = v^ + v^s, s = V°s, v M y^ = o. 



(3.14) 



Now we specialize to three dimensions. The quadratic and gauge fixing part of the action is 
given by 



= i** (□ - f ) (fin + ^±l R + a) fc» 
-^> + C >( D + f)«" + ^[ A ( D + f 



-jV°(° + fl/2) 



A 
R 



1 — 36 r i2i 
-{ D + C+ l} 



2a 



AD + -( □ + - ) I □ + c+ j^j 



A 

^ + 2 



(1 - 36) 2 r i?y , 



(C) 2 " 2g + ^ - W] • (3-15) 



where 



A ee (8a + 3/?)D + 2(3a + /?) fl-<r, 



A 



An 



aR 3a + P 2 
6 18 



Next we choose the parameters a and c such that the action is diagonalized: 

R 2(3a + (3)R/3 - a 



1 - 36 = 2(8a + 3/3)a, c + 



8a + 3/3 



(3.16) 



(3.17) 



We keep the gauge parameter b arbitrary in order to check the gauge dependence of the result. 
Then the action simplifies to 



+ 3(1-36) ^ + \ hAhh > (3 - 18) 
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where we have defined 



A T = (o-|)( /3D + ^ JJ + ,) + 2A, 

A-?2±»e«)(D + f) + (i-»)A, 

A /3-6 i?\ 1-36, 

= A (— D+ 6)-2 A - (319) 

It appears that the choice 6 = 3 makes the kinetic operators simple. However, in this gauge 
the uniform structure of the kinetic terms being quadratic in the D'Alembertian is lost, and the 
method we use does not seem to be suitable to deal with this case. Therefore in the rest of this 
paper we assume that 6 < 3. 
The FP terms (3.13) become 

C = -iiV^AvV,, - ~SA S S) (3.20) 

where 



/ 2(3a + P)R/3 - a R\ / R\ A /3-6 2R\ 

A H° + 8a I W ?)( D + l)' A * = A (— D + t)' <3 ' 21) 

We are now ready to evaluate the beta functions. 



4 Evaluation of the functional traces 

In this section, we set up the calculation of the RG equations for the space of S 3 viewed as 
Euclideanized de Sitter space. 

The quartic structure in the derivatives as opposed to quadratic in [17] can be factorized 
into the product of two quadratic derivatives. Thus we write the kinetic operators in a product 
forms of quadratic derivatives 

^A T = A T1 A T2 = (n-Bl)(n-Bj), (4.1) 

and similarly for Ag, A,,, A^, Ay and As, where 

a + 2Ra ± ^[a + 2i?(3a + /3)/3] 2 - 8/3 A 



B l,2 



2/3 



c a - 2R(3a + /3)/3 ± y/[a + 2i?(5a + 2/3)/3] 2 + 4(36 - l)(8a + 3/3) A 
t>\ - - - - - -' — 



'1,2 



B l,2 



2(8a + 3/3) 
(6 - 3)<r - 2[(36 - 25)q + (6 - 9)/3)]i?/3 



B 



2(6- 3) (8a + 3/3) 

y/[(b - 3)(t - 2{(36 + Tja + (6 + 3)/3}i?/3] 2 + 8(6 - 3)(36 - l)(8a + 3/3) A 
± 2(6-3)(8q + 3/3) ' 

h (6-3)o--2[(36 - 25)a + (6-9)/3)]i?/3 



1,2 



2(6-3)(8a + 3/3) 
^[(6 - 3)<r - 2{(36 + 7)a + (6 + 3)/3}i?/3] 2 - 16(6 - 3)(8a + 3/3)A 
± 2(6- 3) (8a + 3/3) : 
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B 



v 

1,2 



£>1,2 - 



R a + (2a + (3)R/3 
._ 

4i? 



8a + 3/3 

cr - 2(3a + p)R/3 



3(6-3)' 



8a + 3/3 



(4.2) 



For each spin component, we choose the cutoff to be a function of the corresponding operator 
given in (3.19). The gauge fixed inverse propagator is 



o = z 



( \*r 



\ 



4 A 

3(1-36) 



\ 



(4.3) 



The cutoff is chosen as 

/ \n k (A T ) 

n k = z 



\ 



3(1-36) 



K k {A v ) 



(4.4) 



hn k (A h ) J 



We then have 



o + n k = z 



( \Pk{A T ) 



\ 



3(1-36) Pk(A v ) 



(4.5) 



\Pk{A h ) J 



where we have defined the function P k (z) = z + lZ k (z). These formulae are valid for the case 
when the kinetic operators consist of single D'Alembertian, but we should understand that they 
become a sum of those terms when they are products of such operators. We then find 



d t r k = 



l 



Ti{W(A T1 ) + W(A T2 )} + Tr{W(A (1 ) + 1^(A 52 )} + Tr{W(A vl ) + W{A v2 )} 



+Tr{W(A hl ) + W(A h2 )} - Tr{W(A vl ) + W(A V2 )} + Tr{W(A S i) + W(A S2 )} , (4.6) 



where we have defined W(z) 



and t = log k. 



Following [24], we use the optimized cutoff R k (z) = (k 2 — \z\)0(k 2 — \z\). We have d t R k (z) = 
2k 2 6(k 2 - \z\), P k {z) = k 2 for z < k 2 and W(z) = 26(k 2 - \z\). Dividing the numerator and 
denominator by k 2 , they are given by 

d t T k = J> T 0(1 - + J>€*( 1 " I^D + E m ^ " ™ 

n n n 

+ E^(! - \~ X n\) ~ 2X>v*(l - lA^I) - 2^m s e(l - |A*|), (4.7) 



where A^ = \$ /k 2 are the distinct dimensionless eigenvalues of the Euclideanized operator A.;, 
and m$ their multiplicities derived in [17]. The eigenvalues are given by 



a: 
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A nl,2 - 


6 (n 


+ 2n-l) + B* 2 , 


n > 2, 


A nl,2 = 


R, 2 

6 (n 


+ 2n) + B? j2 , n 


> 2 


A nl,2 — 


R, 2 
6 (n 


+ 2n) + Sj 2 , n 


> o 


A nl,2 — 


^/ 2 

6 (n 


+ 2n-l) + B£ 2 , 


n > 1 


A nl,2 — 


R( 2 

6 (n 


+ 2n) + 5^ 2 , n 


> 1 



(4.8) 

where the lower end of the allowed integer n (denoted by no below) for each eigenvalue is also 
shown. Here -B's are those defined in (4.2), but it should be understood that we have rescaled 
the coupling constants by 

G = Gk~\ A = Ak 2 , a = ak~ 2 , p = /%" 2 , (4.9) 

so that they are dimensionless. The multiplicities are 

m T n = 2(n 2 + 2n-3), 
m n = m n = ^(n 2 + 2n), 

ml = m^ = m s n = {n+ I) 2 . (4.10) 
Each sum in (4.7) is evaluated by the Euler-Maclaurin formula 

00 fTlmax 1 f> ft. 

E = I F ^ dn + 2 F{m) - -i F '^ - -i F "'^ + n > ( 4 - n ) 



n=no 



where B n are the Bernoulli numbers (B2 = 1/6 and B4 = 1/30) and 1Z is a remainder. For the 
evaluation of the beta functions, only the first three terms are necessary. Here no are the lower 
ends of the allowed integer n given in Eq. (4.8) and Umax are determined by 

A« = 1. (4.12) 

The evaluation of the sums is done using algebraic manipulation software. As in Ref. [17], 
we restrict ourselves to the parameter region where A and R are of the same order. The reason 
is that the physics is contained in the on-shell effective action though we have to be away from 
on-shell in order to obtain the beta functions. We can be slightly off-shell for this purpose, and 
then we can assume that A and R are of the same order. Therefore we expand the coefficients 
in (4.11) in powers of A and keep at most terms linear in A in A while we keep only the A- 
independent terms in B, and other terms of order AR and R 2 are neglected. Another reason 
why we do not look at the R 2 in this approach is that even if we keep these terms, there is no 
way to tell which are terms corresponding to R 2 and R^, so that we cannot obtain the beta 
functions for a and f3 separately. This is unavoidable as long as we evaluate the functional 
trace for the fixed geometry of sphere. In our approach, the evaluation of the functional trace 
for general background would be difficult because the diagonalization of the quadratic terms of 
gauge-invariant action and gauge fixing term is then more involved and cannot be done explicitly. 

Restricting the terms this way makes the analysis of RG equations enormously simple. In 
our approach we cannot see separately the running of the coupling constants a and /3 of the 
higher derivative terms, which correspond to terms of the order k~ 1 R 2 . Nevertheless we can 
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check the coefficient of these terms. A preliminary study shows that there are fixed points for 
these parameters; we have checked that there are real solutions for a and /3 to the equations 
setting the coefficient to zero. Though this is not sufficient for the existence of the fixed point 
for a and f3, it is a necessary condition. So here we content ourselves with assuming that they 
are fixed at a fixed point in the following discussions. We leave more detailed study of these 
terms for future work. 

The result can be written as 



16vr 



[k A A + kBR + 0{AR)}, 



(4.13) 



where V(S 3 ) = 27r 2 (4) 3 / 2 is the volume of S 3 . The coefficients are found to be 



A 



B 



7T 



{( 



1 + 



3 

<J\ 2 



a 



8a + 3/3 



-2(6-3)Jl + ,-(36-7)6jl 



a 



8a + 3/3 



= A + 4iA, 



3vr(6-3)(8a + 3/3)/3 

+(96 2 - 226 + 3)/3}/3 



a 



(6 - 3) (96a 2 + 28a/3 - 3/3 2 ) I + + {12(26 2 - 56 + 1)5 



P 



8a + 3/3 



= + (36 - 56 - 36 2 )(8a + 3/3)/3 



(4.14) 



Note that there is no term independent of R. This is to be expected for the following reason. 
The result should contain only integer powers of R and in a three-dimensional manifold without 
boundary the volume prefactor is proportional to i?~ 3 / 2 . This implies that the expansion of 
dfT/, contains no i?-independent term. This gives the first nontrivial check of our result. There 
is a term of order R 2 with negative power of k. This is because, as has been noted in [12], there 
is no divergence to the R 2 and R 2 ^ terms in this renormalizable theory. The coefficient of these 
terms is very complicated. 

Evaluating the Euclidean version of the renormalized action (3.1) on the S 3 background, we 

get 



V(S 3 )( 



2A n 



a 



:R 



3« + P 2 



V16ttG 16vrG" 487rG 
Applying the rescaling (4.9) and comparing (4.13) with (4.15), we obtain 



(4.15) 



8vr V g 
a (dtG 
16vr V G 2 



Ad t G 3A 
G 2 + ~G~ 



A 

16tt' 



1 

G 



B 

T&k' 



(4.16) 



Hence we find 



d t G = G(l + aBG), d t A 



-2A + -(A + 2aBA)G. 



(4.17) 



Because we treat the coupling a and /3 as fixed, the RG equations are considerably simplified. 
These equations are exactly of the same form as in pure gravity with cosmological constant and 
in the topological gravity [17]. 
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As a check, we examine if the running of the dimensionless combination 



d t (AG 2 ) = -(A + 6aBA)G 3 , 



(4.18) 



is gauge independent. We find from (4.14) that indeed this is independent of the gauge parameter 
b. 



d t (AG 2 ) = 



3vr 
4 

+ -CJ 

TT 



8a + 3/3 



(8d + 3/3)/3 V P 8a + 3/3 
This gives the second nontrivial check of our results. 



a 



8a + 3/3 



A. (4.19) 



5 Renormalization group flow and fixed points 

We now discuss the properties of the RG equations and fixed points for the theory with a = ±1 
separately. 



5.1 a = +l 

Let us first consider the case when the Einstein term has the usual sign. Our RG equations (4.17) 
are well defined for the range < -1 or ft > 0), and (8a + 3/3 < or 8a + 3/3 > 1). We shall 
be interested in the regions where unitary theories exist; (d ~ — |/3, (3 > 0) or (d > 0, (3 ~ 0). 
So the most interesting region is (d < 0, (3 > 0) and (d > 0, (3 > 0). 

The RE equations have two fixed points. One is the Gaussian fixed point G = A = 0, which 
is seen to be attractive in the A direction and repulsive in the G direction. The other fixed point 
is at 



G* 



1 



A* 



An 



Ai + 6B' 



(5.1) 



where the constants Aq, Ai and B are defined in Eq. (4.14). The constants G and A are dimen- 
sionful parameters, and on general ground we expect the fixed points depend on the gauge. To 
get some idea what values they typically have, we give their values G* ~ 0.133 and A* ~ 0.057 
for d = —1, f3 = 3, b = 0. The flow is shown in Fig. 1 for this choice of d, (3 and b. 

We have checked how much these fixed points are dependent on the gauge parameter b. We 
show how G* changes for —1 < b < 3 in Fig. 2 (a) for d = — 1, /3 = 3 and in (b) for d = 1, (3 = 3. 
These parameters are chosen because we consider that it is natural to choose these of the order 
1. We see that they both give qualitatively similar values and behaviors. In particular the 
absolute value is very small well below b ~ 3, so that this is certainly within the perturbative 
domain. Therefore, though our result is based on a onedoop calculation, we can expect that it 
correctly represents the features of the Wilsonian RG flow. 

It is remarkable that the 6-dependence disappears from the fixed point of the cosmological 
constant. As we saw in the preceding section, only A\ and B depend on b, but the combination 
A\ + 6B is precisely that appearing in the beta function of the dimensionless coupling (4.18), 
and hence that dependence drops out. We find that for the typical choice of d and (3, the 
cosmological constant is positive and small; 0.057 for a = — 1 and (3 = 3, and 0.045 for d = 1 
and (3 = 3. The precise values may not be so important, but the fact that they are generally 
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rather small may have some significance. Fig. 3 shows how the values of G and A change for 
< a < 3 and < j3 < 3. This figure shows that the values are indeed small for a wide range 
of a and /3 except near 8a + 3/3 = and (3 = 0, where these quantities are singular. 

5.2 a = -l 

Let us now turn to the opposite sign for the Einstein term. Our RG equations (4.17) are well 
defined for the range < or (3 > 1), and (8a + 3,5 < -1 or 8a + 3/3 > 0). We shall be 
interested in the regions where unitary theories exist; (a « — 1/3,/3 > 0) or (a > 0, j3 « 0). So 
the most interesting region is (a < 0, (3 > 1) and (a > 0,/3 > 0). 

Here again there is the Gaussian fixed point which is attractive in the A direction and 
repulsive in the G direction. Another nontrivial fixed point in this case is given by 
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(a) (b) 
Figure 3: The fixed point values of (a) G* and (b) A* for b = as a function of a and j3. 




Let us see again the concrete values and the gauge dependence of the fixed point of the grav- 
itational constant. For a = —l,/3 = 3, b = 0, we find G> ~ —0.161 and A* ~ 0.106. The flow 
is shown in Fig. 4. We show how changes for — 1 < b < 3 in Fig. 5 (a) for a = — 1,/3 = 3 
and in (b) for a = l,/3 = 3. We see that again they both give qualitatively similar values and 
behaviors. In particular the absolute value is again very small. 

We note that the sign of the Newton constant changes from a = +1 case and it takes negative 
value for typical values of a and f3. Recall that in this case we took the negative sign for the 
Einstein term, but the fixed point of the gravitational constant takes negative values, resulting 
in positive Einstein term. Namely the Einstein term has positive coefficient at the fixed point 
even if we start with negative sign. 

It is remarkable that the fixed point of the cosmological constant again does not depend on 
the gauge. We find that its values are 0.106 for a = —1, (3 = 3 and 0.038 for a = 1, j3 = 3, again 
very small positive numbers. Fig. 6 shows how the values of G and A changes for —3 < a < 3 
and < p < 3. 
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(a) (b) 
Figure 6: The fixed point values of (a) and (b) A* for b = as a function of a and f3. 



6 Discussions and conclusions 

In this paper we have studied the quantum effects of three-dimensional higher derivative gravity 
which has attracted much attention recently. In particular we have obtained the Wilsonian 
RG equations to study their fixed points. Though we did not attempt to derive those for the 
coefficients of the higher curvature terms because we cannot derive those for a and j3 separately 
and also they become awfully complicated, we have some evidence that there is some fixed point 
for a and f3. There are several conclusions we can draw even within this restriction. Assuming 
that a and j3 have fixed points, we have found that there are ultraviolet fixed points for the 
gravitational and cosmological constants, one the usual Gaussian and the other nontrivial one 
for both signs of the Einstein term. This shows that this theory is asymptotically safe. 

We have also found that the fixed point value of the gravitational constant is generically 
small, so that our one-loop discussion may well be justified. What is interesting is that if we 
take the negative sign for the Einstein term in the bare action, the gravitational constant has 
only a fixed point of negative values or zero, in the former case resulting in positive Einstein 
term. This is an interesting result which has not been explored in other approaches. It would 
be extremely interesting whether this happens also in four-dimensional gravity. 

We have also shown that RG equations are singular for the parameters, 8d + 3/3 = or (3 = 0, 
corresponding to the new massive gravity or a special case of f(R) gravity. Indeed, we see from 
Figs. 3 and 6 that and A* diverge there. If we look at the diagonalized action (3.18), we find 
that only the kinetic term for the transverse traceless mode is quadratic in the D'Alembertian 
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but the rest are all linear, as can be seen from (3.19), and the limit is singular. This means that 
the new massive gravity and f(R) gravity would not correspond to a fixed point in this class 
of theories. We would like to note that this result is one of the conclusions that are obtained 
within the approximation that the coefficients of the higher curvature terms a and f3 are not 
subject to the flow, so strictly speaking there may be a possibility that the picture might be 
different in a larger space when these parameters are allowed to run. 

Thus, though we did not attempt to evaluate the beta functions for a and $, it is quite 
interesting to examine whether the picture may change if we include those parameters in the 
RG analysis. For this purpose, we have to find a way to compute these for general backgrounds, 
not just for sphere, because we have to tell which is the contribution to R 2 and which is to R^ v . 
Probably the heat kernel method may be useful in such a study. We hope to return to this 
problem in the future. 
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